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ABSTRACT 

We investigate the effects of minor mergers between an SO galaxy and a gas-rich satellite galaxy, by means of N- 
body/smoothed particle hydrodynamics (SPH) simulations. The satellite galaxy is initially on a nearly parabolic orbit 
and undergoes several periapsis passages before being completely stripped. In most simulations, a portion of the stripped 
gas forms a warm dense gas ring in the SO galaxy, with a radius of ~ 6 — 13 kpc and a mass of ~ 10' Mg. The ring 
is generally short-lived ( 3 Gyr) if it forms from prograde encounters, while it can live for more than 6 Gyr if it is 

born from counter-rotating or non-coplanar interactions. The gas ring keeps memory of the initial orbit of the satellite 
galaxy: it is corotating (counter-rotating) with the stars of the disc of the SO galaxy, if it originates from prograde 
(retrograde) satellite orbits. Furthermore, the ring is coplanar with the disc of the SO galaxy only if the satellite’s orbit 
was coplanar, while it lies on a plane that is inclined with respect to the disc of the SO galaxy by the same inclination 
angle as the orbital plane of the satellite galaxy. The fact that we form polar rings as long-lived and as massive as 
co-planar rings suggests that rings can form in SO galaxies even without strong bar resonances. Star formation up to 
0.01 M@ yr -1 occurs for > 6 Gyr in the central parts of the SO galaxy as a consequence of the interaction. We discuss 
the implications of our simulations for the rejuvenation of SO galaxies in the local Universe. 

Key words. Galaxies: interactions - Methods: numerical - Galaxies: kinematics and dynamics - Galaxies: ISM 
galaxies: star formation 


1. Introduction 

A large fraction of SO galaxies (~ 20 — 25%, 

lLaurikaineu et al.ll201.3D show signatures of rings (i.e. closed 
structures composed of stars, gas and dust) and pseudo¬ 
rings (incomplete versi ons o f ringSj_s omet imes formed 

by arm-like structures, IComeron et al.l 12014 ). Far ultra¬ 
violet (FUV) observations frequently evidence rings and 
pseudo-rings of gas and/or young stars in early -type galax¬ 
ies (ETGs), and especially in SO galaxies lThilker_et - ^l 
l2007t ISalim fe Richll2010l : ISalim et~al1l2012D . ISalim fe Rich 
( 20101) analyze 22 ETGs with an extended ultraviolet (UV) 
morphology in GALEX data, and find rings in 15 of them. 
ISalim etahl (120121) investigate a sample of ETGs at z < 0.1, 
27 of which have been observed with the Hubble Space 
Telescope ( HST) in FUV. These ETGs are UV structured, 
showing ring and/or arm-like structures, while in the opti¬ 
cal they have an old stellar disc, i.e. they are SO galaxies. 
The UV emission is consistent with a low-level galaxy-scale 
star formation (SF) of ~ 0.5 M 0 yr -1 . 

Rings in SO galaxies are found at different distances 
from the galaxy n ucleus, often connec t ed to the presence 
of ba r/s (see e.g. lAguerri et all 120091 : lLanrikainen et all 
120131) . Such kind of rings are often classified, depending on 
their distance from the nucleus, as (i) nuclear rings (sur¬ 
rounding the nucleus of the galaxy), (ii) inner rings (gen¬ 
erally surrounding the bar, if present), and (iii) outer rings 
(more external than the former ones). 


The origin of the rings is still debated. The fact 
that many rings are associated with bars (or other non- 
axisymmetric structures) has led to the idea that the for¬ 
mation of rings is connected with resonances induced by the 
rotation of the bar (i.e. either an inner Linblad’s resonance, 
or an o uter Linblad’s resonance, or a ultraharmonic res - 
onance, iBinnev fc Tremainel 1 1 987t iButa fc Gombesl[l99fih . 
However, a re sonance can n ot be invoked for the outer po¬ 
lar rings (e.g. iMarino et alll2009h . as well as for the mos t 
asymmetric and off-centre rings (e.g. iThilker et al.l 1201011 . 
In addition, bars shoul d be more co mmon than observed 
in SO galaxies (~ 30%. lAguerri et, alll2000fl to explain the 
l arge fraction of ringed SO galaxies found in ISalim fc Rlchl 
( 120101 ) . 

Alternative scenarios to a bar-driven resonance include 
gas-rich minor or major mergers, and accretion from fil- 
aments/clouds of intergalactic medium. Mergers and gas 
accretion from filaments are the only models able to ex- 


plain the most massive polar rings (Bekki 19981: 

odice et al.. 

2002, 20061 Moiseev et al. 20 lit 

Spavone et al. 

201(1 2011: 

Soavone & Iodice 2013|). Evidences of gas accretion either 


from filaments or from interactions with gas-rich satel¬ 
lite gala xies have been found in a large number of ga lax¬ 
ies (e.g. Ide Blok et al.l 12014 ISerra fe Oosterlool 120101 : see 
ISancisi et al.l 120081 for a review). Similarly, the analysis 
of a large sample of d isc galaxies in the Sloan Digital 
Sky Survey reported by iKavirail (|2014af ) suggests that mi¬ 
nor mergers played a major role in driving SF in galaxy 
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Table 1 . Initial conditions of the TV—body simulations: 
masses and scale lengths. 


Model galaxy properties 

Primary 

Secondary 

Mdm (10 11 Mq) 

7.0 

0.3 

M, (10 10 M 0 ) 

7.0 

0.2 

/b/a 

0.25 

0.25 

M g (10 8 M 0 ) 

0 

1.38 

R s (kpc) 

6.0 

3.0 

Rd (kpc) 

3.7 

3.0 

hd (kpc) 

0.37 

0.30 

r B (kpc) 

0.6 

0.6 


Mdm is the DM mass; M* is the total stellar mass of the 
galaxy (including both bulge and disc); fb/d is the 
bulge-to-disc mass ratio; Mg is the total gas mass. The 

primary has no gas, while the gas of the seconda ry is _ 

distributed according to an exponential disc (iHernauistll 19931 . 
with the same parameters (scale length and height) as the 
stellar disc. R B is the NFW scale ra dius R s = -R 200 /c, wh ere 
R 200 is the virial radius of the halo (iNavarro et all . 119961 1 and c 
the concentration (here we assume c = 12 for both galaxies). 

Rd and hd are th e disc scale length and height, respectively 
(lHernouistlll993 h. while r B is the bulge scale length. 


discs (see also iKavirail 120101 : iKavirai et al.l l2012t 120131 : 
iKavirail l2014bl) . Numerical simulations also indicate the 
importance of mi nor mergers for the accretion of fres h 
gas in E TGs (e.g. Kavirai et al.ir2009HPeirani et a.]1l2010t) . 
Finally, IMarino et al. (|201lD investigate five SO galaxies 
with outer rings in UV (NGC 1533, NGC 2962, NGC 
2974, NGC 4245, and NGC 5636). While all the consid¬ 
ered galaxies are barred and the rings might be related to 
bar-induced resonances, there are several indications that 
NGC 1533, NGC 2962 and NGC 2974 ha v e acc reted gas 
from satellite galaxies. Thus, iMarino et al.l (120111) hypoth¬ 
esize that the outer rings in their sample are an effect 
of (bar-driven) secular evolution, but at the same time 
this process was aided by fresh gas accreted from gas-rich 
satellites. Using N-body/smoothed-particle hydrodynamics 
(SP H) simulations with cliemo-photometric implementa¬ 
tion, iMazzei et al.l (|2014l ) showed that the pseudo-ring in 
NGC 1533 and its star forming properties may be a tran¬ 
sient phase during a major merger episode involving two 
halos with mass ratio 2:1. 

In this paper, we study the formation of gas rings in SO 
galaxies, by means of N-body/SPH simulations of minor 
mergers with small gas-rich companions. In Section 2, we 
describe the adopted numerical techniques. In Section 3 
we present our results, with particular attention for the 
formation of co-planar and polar gas rings. In Section 4 
we discuss the main implications of our results, while in 
Section 5 we summarize our conclusions. 


2. Method: AT—body simulations 

We simulate interactions between a primary SO galaxy 
and a secondary gas-rich galaxy. The initial conditions for 
both the primary galaxy and the secondary galaxy in the 
TV—body model are generated by using an upgraded version 
of the code d escribed i n Wi drow, Py m fe D ubinski J(2008; 
see a lso iKuiiken & Dubinskilll995l and IWidrow fc Dubinsldl 
l2005f h The code generates self-consistent disc-bulge-halo 
galaxy models, derived from explicit distribution functions 


for each component, which are very close to equilibrium. 
In particular, the halo is modelled with a Navarro, Frenk 
& White (1 996, NFW ) prof ile. We use an exponential 
disc model (iHernauistl . ! 199 31 ). while the bulge is spher- 
i cal and comes from a generalization of th e Sersic law 
dPrugniel &; Simier]ll997nWidrow et al.ll2008l) . 


Both the primary and the secondary galaxy have a stel¬ 
lar bulge and a stellar disc. The giant SO galaxy has no 
gas, whereas the secondary galaxy has an initial gas mass 
of 1.38 x 10 8 M 0 , distributed according to an exponential 
disc. Therefore, the initial configuration of the secondary 
galaxy is consistent with a low-mass gas-rich disc galaxy. 
The assumption that the primary has no gas is quite un¬ 
realistic, but allows us to easily distinguish between the 
triggering of SF in the pre-existing gaseous disc of the SO 
galaxy and the effect of fresh gas accretion from the com¬ 
panion galaxy. 

The total mass of the secondary is ~ 1/20 of the mass 
of the primary, classifying the outcome of the interaction 
as a minor merger. The masses of the various components 
and the scale lengths of the simulated galaxies are listed 
in Table 1. Table 2 shows the orbital properties (impact 
parameter, relative velocity, orientation angles and total 
energy) of the six runs. The adopted orbits are nearly 
parabolic (Table 2), in a greement with pred i ctions from 
cosmological simulations dKochfar fe Burkertl 120061) . The 


main implication of such nearly parabolic orbits is that the 
stripping of the satellite galaxy proceeds slowly: the merger 
is not yet complete in most simulations at t = 11 Gyr since 
the first periapsis passage. 

In all the simulations, the particle mass in the primary 
galaxy is 2.5 x 10 5 M 0 and 5 x 10 4 M 0 for dark matter (DM) 
and stars, respectively. The particle mass in the secondary 
galaxy is 2.5 x 10 4 M 0 for DM and 5 x 10 3 M 0 for both 
stars and gas. The softening length is e = 100 pc. 

We simulate the evolution of the models w i th th e 
N— body/SPH tree code gasoline (|Wadslev et al.l . I2004D . 
Radiative cooling, SF and supernova blastwave feedback 
are enabled, as described in Stinson et al. (2006, 2009, 
see also lKatz fe R.ixlfl992D . The adopted parameters for SF 
and feedback are the same as used in recent cosmological 
simulations capable o f forming realistic galaxies in a wide 
range of masses (e.g., iGovernato et aLlfeolQl : iGuedes et al.l 
[2011 ), and in recent simulations of galaxy-galaxy c ollisions 
(jManelli fo Maveill20 1 2l : iManelli et, al.l 1201 2 . l2013alfbh . 

In runs A and B the orbit of the satellite galaxy is co- 
planar with the disc of the SO galaxy. The only difference 
between run A and B is that the orbit of the satellite galaxy 
is retrograde and prograde (i.e. the orbital angular momen¬ 
tum of the satellite is counter-aligned and aligned with re¬ 
spect to the spin of the primary galaxy) in run A and B, 
respectively (Table 2). Run A, C and D differ for the incli¬ 
nation angle between the relative velocity vector v re i and 
the symmetry axis of the primary disc (0, see Table 2): 
9 = —7r/2, —7r/4 and 0 for run A, C and D, respectively. 
This means that the orbit of the satellite galaxy is inclined 
by 45° and 90° with respect to the disc of the primary 
galaxy in run C and D, respectively. 

In runs E (retrograde) and F (prograde), the impact 
parameter (30 kpc) is a factor of 3 larger than in run A. 
In all simulations, the primary galaxy develops a strong 
bar ~ 200 Myr after the beginning of the simulation. Our 
galaxy models are marginally stable against the formation 
of a bar: evolving the galaxy in isolation, the bar does not 
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Table 2. Initial conditions of the IV—body simulations: orbital properties. 


Run 

b 

(kpc) 

^rel 

(km s -1 ) 

e, <p, ip 

(rad) 

E a 

(10 4 km 2 s" 2 ) 

L s 

(10 3 km s 1 kpc) 

e 

Orbit spin 

A 

10.0 

200 

—7t/2, 0, 0 

0.38 

2.0 

1.003 

retrograde 

B 

10.0 

200 

7t/2, 0, 0 

0.38 

2.0 

1.003 

prograde 

C 

10.0 

200 

—tt/4, 0, 0 

0.38 

2.0 

1.003 

- 

D 

10.0 

200 

0 , 0 , 0 

0.38 

2.0 

1.003 

- 

E 

30.0 

200 

-tt/2, 0 , 0 

0.38 

6.0 

1.03 

retrograde 

F 

30.0 

200 

tt/2, 0 , 0 

0.38 

6.0 

1.03 

prograde 


b and v re i are the impact parameter and the relative velocity bet ween the centres of ma ss (CMs) of the two galaxies at the initial 
distance, respectively. For the definition of 6, <p, ip, see figure 1 of [Hut fe Bahcall (1 1 9831 1 . In particular, 6 is the angle between the 
relative velocity vector v re i and the symmetry axis of the primary disc, cp describes the orientation of v re i projected in the plane 
of the primary disc and ip describes the orientation of the initial distance vector D (between the CMs of the two galaxies) in the 
plane perpendicular to v re i. 

E a is the specific orbital energy, i.e. the total energy divided by the reduced mass fi = mi m 2 /(mi + m 2 ), where mi and m 2 are 
the mass of the primary and of the secondary galaxy, respectively. E s = — G M/D + v 2 el /2, where M = mi + m 2 is the total 
mass of the two galaxies, G is the gravitational constant and D the initial distance between the CMs. 

I/ a is the modulus of the specific orbital angular momentum, i.e. the angular momentum divided by the reduced mass, 
e is the eccentricity (e = [1 + 2 E s L a /(G M) 2 ] 1 ^ 2 ). 

A orbit is classified as prograde/retrograde depending on the alignment/counter-alignment of the orbital angular momentum of 
the secondary galaxy with respect to the spin of the primary galaxy. 



t [Gyr] 


Fig. 1. Radial distance of the satellite galaxy from the cen¬ 
tre of the primary galaxy for all simulations, as a function 
of time. Time t = 0 corresponds to the first periapsis pas¬ 
sage. Long-short dashed black line: run A; solid red line: 
run B; dotted ochre line: run C; dashed green line: run D; 
long-dash-dotted blue line: run E; short-dash-dotted ma¬ 
genta line: run F. 

appear. The perturbation induced by the interaction with 
the intruder satellite galaxy plays an important role in the 
formation of the bar. 

3. Results 

Fig- CD shows the radial distance of the satellite galaxy from 
the centre of the primary galaxy for all simulations, as a 
function of time (in this Figure and in the following, time 
t = 0 corresponds to the first periapsis passage). From this 



t [Gyr] 


Fig. 2. Mass of the gas that lies in the inner 15 kpc of 
the SO galaxy (M gi ; n ) for all simulations, as a function of 
time. Time t = 0 corresponds to the first periapsis passage. 
Long-short dashed black line and filled black circles: run A; 
solid red line and open red circles: run B; dotted ochre line 
and open ochre squares: run C; dashed green line and filled 
green squares: run D; long-dash-dotted blue line and blue 
asterisks: run E; short-dash-dotted magenta line and filled 
magenta triangles: run F. 

figure, we can derive information about the periapsis pas¬ 
sages. The impact parameter is extremely important to de¬ 
termine the amount of tidal disruption of the satellite at 
the first periapsis passage, and the subsequent orbital de¬ 
cay. In fact, the orbital decay is slower in the two runs with 
the largest impact parameter: in 10 Gyr, the satellite un¬ 
dergoes only 2 and 3 periapsis passages in runs E and F, 
respectively. In the other runs, the satellite undergoes be- 
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t [Gyr] 


Fig. 3. Mass of the gas ring (M r i ng ) for all simulations, as 
a function of time. Time t = 0 corresponds to the first peri- 
apsis passage. Long-short dashed black line and filled black 
circles: run A; solid red line and open red circles: run B; dot¬ 
ted ochre line and open ochre squares: run C; dashed green 
line and filled green squares: run D; short-dash-dotted ma¬ 
genta line and filled magenta triangles: run F. 


tween 4 and 8 periapsis passages in 10 Gyr. Generally, the 
orbital decay is faster for prograde satellites (runs B and 
F) than for their retrograde analogs (runs A and E): tidal 
disruption occurs faster in prograde encounters, because 
of the lower relative velocity between the stars of the pri- 
mary galaxy and th e gas a nd stars of the secondary galaxy 
(jToomre fe Toomrei . Il972f) . 

3.1. The gas rings 

At each periapsis passage, gas and stars are stripped away 
from the satellite galaxy, and remain bound to the primary 
galaxy. In some runs, a portion of the stripped gas settles 
into a warm gas ring (T ~ 10 4 K), with an average radius 
of ~ 6 — 13 kpc. Figure [2] shows the total mass of stripped 
gas that resides in the inner 15 kpc of the primary galaxy 
(M g) in), as a function of time. By definition, M gi - m includes 
all the gas that lies in the inner 15 kpc of the SO galaxy, 
including diffuse hot gas, the gas that is funnelled in the 
nucleus of the galaxy, and the gas that is confined into the 
ring. M gi j n becomes non-negligible (^,2x 10 6 M 0 ) soon 
after the first periapsis passage (t ~ 0.5 — 1.0 Gyr) in all 
simulations apart from run E. In run E (retrograde and 
with a large impact parameter), gas stripping starts only 
in the last Gyrs. In all simulations apart from run E, M g: j n 
reaches a maximum value of 1 — 2 x 10 7 M 0 (i.e. ~ 1/10 of 
the total gas mass of the satellite galaxy). From Figure [2j 
we also notice that the peak value of M gt j n is reached earlier 
in the prograde runs (run B and F) than in the retrograde 
runs (runs A and E) and in the non-coplanar runs (runs C 
and D). This confirms that gas stripping is more efficient 
in the prograde runs. The fact that retrograde runs and 
non-coplanar runs behave in a similar way is remarkable. 


Figure [3] shows the total mass of stripped gas that is 
confined in a gas ring (M r i ng ), as a function of time. During 
the ring’s life, the value of M r i ng spans from 20% to 90% of 
M gi i n , depending on the run. From FigurejSJ it is apparent 
that the formation of a ring, its survival, and the amount 
of mass that is locked into the ring strongly depend on the 
characteristics of each run. In some runs, the ring never 
forms (run E), or is relatively short-lived (1 — 3 Gyr, runs B, 
C and F), while in other runs the ring grows steadily and 
is long-lived (> 7 Gyr, runs A and D). 

In run A (Fig. [4]) a co-planar ring of gas forms in the 
disc of the SO galaxy at t ~ 3 Gyr, i.e. approximately after 
the second periapsis passage. The ring has a radius of ~ 
8 — 10 kpc, which means that it forms just at the end of 
the stellar bar, but is counter-rotating with respect to the 
disc, because it keeps memory of the retrograde orbit of 
the satellite. The ring in run A is extremely long-lived: it 
survives for the entire simulation (Figure [3]). In contrast, a 
gas ring forms earlier (at t ~ 2 Gyr) in run B (the prograde 
analogous of run A), is co-rotating with the stellar disc, and 
survives only for ~ 1 Gyr (Figure [3]). The main reason for 
the early disappearance of the ring in run B is that most 
of the gas is channelled inwards by the bar, and ends up in 
the nucleus. 

The fact that the radius of the ring coincides almost 
with the size of the bar in runs A and B might lead to the 
conclusion that the ring is the consequence of a bar-driven 
resonance. On the other hand, rings form even in runs C and 
D (Figure[5]), in which the orbit of the satellite is inclined by 
45 and 90° with respect to the disc of the primary galaxy, 
respectively. Actually, it seems that the ring preserves the 
inclination of the orbit of the satellite with respect to the 
primary galaxy, as the gas ring of run C and D forms with 
an inclination of ~ 45 and ~ 90° with respect to the disc of 
the primary galaxy. In particular, the gas ring of run D is 
a true polar ring, and cannot be produced by a bar-driven 
resonance. The ring of run D grows steadily in mass and 
survives till the end of the simulation. The radius of the 
ring in both run C and D is ~ 12 — 13 kpc, slightly larger 
than in the other runs. 

The fact that the ring has approximately the same ori¬ 
entation as the initial satellite orbit is consistent with our 
expectations: no dynamical process (e.g. bar-driven reso¬ 
nance, differential precession) is sufficiently efficient to per¬ 
turb the initial orbit of gas and stellar streams within a 
few Gyr, if the satellite’s orbit has a large incli nation . 
Even shocks and virial heating (e.g. Ilveres et~aD l2005h 
are partially inefficient for our systems, giv en the lack 
of gas in the halo of the SO galaxy. Finally, iTheis et al.l 
( 2006 1 find that spiral arm-like perturbations might de¬ 
velop in a polar ring surrounding a SO galaxy, but these 
perturbations are confined in the initial plane of the polar 
ring. We warn that our results are valid only for spher¬ 
ical halos. A triaxial halo is expected to exert torques, 


Suarkd 1984 

David et al. 

1984 Steiman-Cameron et al. 

1984 

: Nelson & Tremaine 

199 * 3 : Dubinski & Chakrabartv 


200 91. Such effects might perturb a polar gas ring signif¬ 


icantly. 

Furthermore, our findings are consistent with previ¬ 
ous studies which investigate the formation of polar rings 
from the disruption of gas-rich sat ellite galaxies (e.g. 
iRichter et al.lflOOltlMazzei et al.ll2014H . In particular, there 
is consensus that the initial plane of the ring is the or- 
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Fig. 4. Projected density of gas (logarithmic colour-coded map) and stars (isocontours) of the primary galaxy in run A. 
From top left to bottom right: 1, 2, 2.35, 3, 5, and 7 Gyr after the first periapsis passage. Each panel measures 79 x 67 
kpc. In the top right-hand panel, the white circle marks the position of the nucleus of the satellite galaxy, ~ 50 Myr after 
the second periapsis passage. 


bital plane of the companion if this is prompt ly destroyed 
(e.g. iKatz fe Rbd[l992t iThakar fe Rvdenlll996lf . The initial 
ring plane might then evolve to a preferred plane (either 
the equatorial or a polar plane) on a timescale longer than 
the orbital or azimuthal smearing time (i .e. sever al Gyr, 
e.g. iRix fc Kat7l[T99i1: lG ox fc Snarkelll99fl see also [Struck] 
I1999L and references therei n). Based on a sample of 78 po¬ 
lar ring galaxy candidates, ISmirnova fe Moiseevl (1201311 re¬ 
cently proposed that the DM halo and the stellar bar are 
crucial to stabilize outer and inner polar rings, respectively. 
This is fairly consistent with our results, which suggest that 
the formation of the ring is not due to bar-driven resonance, 
but the bar and the DM halo contribute to stabilize the 
ring, and make it long-lived. 

In both runs E (retrograde) and F (prograde) the orbit 
of the satellite is co-planar with the disc of the primary 
galaxy. The main difference with runs A and B is that the 
periapsis passage is 3 times larger in runs E and F. Only in 
run E does no significant ring form over the duration of the 
simulation (t ~ 10 Gyr), because the periapsis distance of 
the satellite is too large and the relative velocity too high 
to produce significant stripping. 

In run F, which is the prograde analogous to run E, 
a ring forms already at t ~ 3 Gyr since the first periapsis 
passage (Figure[3|. This confirms that gas stripping is much 
more efficient if the satellite orbit is prograde. Despite the 
larger periapsis with respect to the other runs, the ring of 
run F forms again at the end of the stellar bar, with a radius 
of ~ 10 — 12 kpc. The ring of run F disappears quite fast 
(t ~ 6 Gyr). This result is analogous to what we already 


discussed for run B: in the prograde runs the ring appears 
earlier than in the retrograde ones, but disappears earlier. 

3.2. The long-lived gaseous and stellar halo 

A considerable fraction of the initial stellar and gaseous 
mass of the satellite galaxy is scattered into the halo of the 
SO galaxy, as a consequence of the encounter. Fig. [6] shows 
the total mass of stars (M* )OUt ) and gas (M giOUt ) that was 
initially bound to the satellite galaxy and that, at time t, 
lies at a distance > 50 kpc both from the centre of the 
SO galaxy and from the centre of the satellite galaxy. In 
our definition of M*. out and M giOUt , we exclude all stars 
and gaseous particles that he at a distance < 50 kpc from 
both the nucleus of the SO galaxy and the satellite galaxy, 
because we want to estimate the mass of the diffuse stel¬ 
lar and gaseous halo, which forms as an effect of the tidal 
encounter (a distance of 50 kpc is a conservative limit, to 
exclude matter that is still in the disc of the SO galaxy or 
that can be re-accreted by the satellite galaxy). 

From Fig.[6]it is apparent that Af* iOUt and M gjOUt grow 
very fast after the first periapsis passage (the only excep¬ 
tion being run E). After the first Gyr, M* jOUt and M giOUt 
remain almost constant, or keep growing slowly, till the end 
of the simulation. Thus, the first periapsis passage plays a 
major role in the formation of the stellar and gaseous halo. 
Furthermore, the halo is extremely long-lived. At the end 
of the simulation (10 Gyr after the first periapsis passage), 
~ 70 — 95% of the stars initially bound to the satellite 
galaxy lie in the stellar halo of the SO galaxy (including the 
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Fig. 5. Projected density of gas (logarithmic colour-coded 
map) and stars (isocontours) of the primary galaxy 4 Gyr 
after the first periapsis passage. Top panel: run C. Bottom 
panel: run D. Each panel measures 79 x 67 kpc. 


case of run E). Similarly, ~ 60 — 80% of the gas initially 
bound to the satellite galaxy ends up in the gaseous halo of 
the SO galaxy at t = 10 Gyr. The only exception is run E, 
where M g)OUt is only ~ 30% of the initial gas mass at the 
end of the simulation. 

3.3. The stellar shells 

The main feature of the minor mergers in the stellar com¬ 
ponent of the SO galaxy is the formation of shells (FigurcQ. 
Shells are ubiquitous in the stellar component of our simu¬ 
lations: they start forming already during the first apoapsis 
passage, and last for the entire simulation in all runs. Our 
findin g is consisten t with the results of previous studies (e.g. 
iPeirani et al.ll2010h . 

In all simulations, the mass confined in the stellar shells 
is ^ 0.5 M* jOUt , i.e. it is more than half of the total mass 
of the stellar halo discussed in Section EQ1 (Fig. l6l) . On the 
other hand, the surface brightness of most shells is quite 



0123456789 10 
t [Gyr] 


Fig. 6. Top (Bottom) panel: mass of stars (gas) initially 
bound to the satellite galaxy that lie(s) at a distance > 
50 kpc both from the centre of the SO galaxy and from 
the centre of the satellite galaxy, as a function of time. 
M*, 

out (M g , out) can be considered as the mass of the stellar 
(gaseous) halo of the SO galaxy. Time t = 0 corresponds to 
the first periapsis passage. Long-short dashed black line and 
filled black circles: run A; solid red line and open red circles: 
run B; dotted ochre line and open ochre squares: run C; 
dashed green line and filled green squares: run D; long- 
dash-dotted blue line and blue asterisks: run E; short-dash- 
dotted magenta line and filled magenta triangles: run F. 

low. For example, the total mass of stars in the shells in 
run A at 5 Gyr is ~ 7 x 10 8 M 0 , i.e. it is nearly 35% of 
the initial stellar mass of the satellite galaxy, but this mass 
is spread over ~ 400 x 400 kpc. The maximum density in 
the strongest shells visible in Fig.[7]is ~ 1.4 — 2.6 x 10 4 Mg 
kpc -2 , much lower than the density in the outer parts of 
the SO disc (~ 10 7 M 0 kpc -2 ). At 7 Gyr, the total mass of 
stars in the shells in run A rises up to ~ 9 x 10 8 M 0 , i.e. it is 
nearly 45% of the initial stellar mass of the satellite galaxy, 
but the density in the single shells and the ratio between 
density of the stellar disc and density of the shells do not 
change. Most shells in our simulations lie at much larger 
radii than the disc scale-length of the SO galaxy, and are 
composed of matter removed from the satellite galaxy. We 
do not find evidence of shells com posed of matter fronmthe 
SO galaxy (as proposed bv I Wallin fc Struck-Marceilll988f) . 

3.4. The star formation rate 

Our simulations include recipes for stochastic SF, as de¬ 
scribed in Section [2] Because of our choice of not including 
gas in the initial conditions of the SO galaxy, only the gas of 
the satellite galaxy can contribute to SF. This allows us to 
study SF from stripped gas in detail. Fig. [8] shows the SF 
rate as a function of time in our simulations. In most simu¬ 
lations, SF starts at the first periapsis passage (t = 0), and 
occurs mainly in the centre of the satellite galaxy, which 
has been strongly perturb ed by the close interaction (see 
iManelli et al .112012L l2013al for more details). As soon as gas 
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Fig. 7. Projected density of stars (logarithmic gray-scale 
map) 5 Gyr since the first periapsis passage in run A. The 
box measures 400 x 400 kpc. The stellar disc of the SO (seen 
face-on) is saturated. 



t [Gyr] 


Fig. 8. Star formation rate in all simulations, as a func¬ 
tion of time. Time t = 0 corresponds to the first periapsis 
passage. Long-short dashed black line: run A; solid red line: 
run B; dotted ochre line: run C; dashed green line: run D; 
long-dash-dotted blue line: run E; short-dash-dotted ma¬ 
genta line: run F. 


is stripped and accretes onto the SO galaxy, SF takes place 
mainly in the central parts of the SO galaxy. Thus, the first 
SF burst in Fig. [5] occurs in the satellite galaxy during the 
first periapsis passage, while the second, longer SF episode 
occurs at the centre of the SO galaxy, and is powered by 
stripped gas. We notice that the SF episode in the satellite 
galaxy and the SF episode in the centre of the SO galaxy 


have a very different duration: they last for < 1 Gyr and 
for ~ 8 Gyr, respectively. 

Interestingly, the SF rate is higher in runs B, F (i.e. 
the prograde runs) and C. These are also the three runs 
where the ring is more short-lived (apart from run E, in 
which the ring does not form). This suggests that the ring 
is suppressed when SF forms stars and stars explode as 
supernovae heating the surrounding gas. On the other hand, 
in our simulations most SF occurs in the nucleus of the SO 
galaxy, rather than in the ring. This is likely an i ntrinsic 
limit of the stochas tic recipes for SF (see e.g. (Mapelli et all 
I2008L I2012L l2013alf : these cannot resolve SF in relatively 
diffuse warm gas regions (the ring), while enhance SF in 
the density peaks (the nucleus). 

For most simulated galaxies, the SF rate is ~ 10 -2 — 
10 -3 Mq yr -1 , almost constant for several Gyr. The run 
where gas stripping is less efficient (run E) is also the run 
with the lowest SF rate (~ 10 -4 M 0 yr -1 , one order of mag¬ 
nitude lower than in the other runs). Our simulated SF rate 
is quite lower tha n the SF rate o f rejuvenated SO galaxies 
as estimated bv I Salim et al.l (l2012f) (~ 0.5 M 0 yr -1 ). This 
difference can arise from the fact that our simulated SO 
galaxies are completely gas-free before the merger, while 
real SO galaxies have a gas reservoir, in which SF can be 
triggered during the merger. 


4. Discussion 


Gas rings form in most of our simulations. Thus, minor 
mergers with gas-rich satellites appear to be a very promis¬ 
ing scenario for the formation of the gas rings observed in SO 
galaxies. But how many of these SO rings can be explained 
by minor mergers? 

Cosmological simula tions co mbined with semi-analytic 
models (e.g. iBertone fc Conselicel l2009i f indicate that the 
minor-merger fraction at redshift z = 0 (i.e. the fraction of 
galaxies that undergo a minor merger at z = 0) is / mm ~ 
0.1 — 0.3, for galaxies with stellar mass M* ~ 0.3 — lx 10 11 
M 0 (a reasonable mass range for SO galaxies). Since the 
merger-induced SF episodes in our simulations are long- 
lived ( > 6 Gyr, Fig.®, we expect that all SO galaxies that 
undergo a minor merger at z = 0 show signs of rejuvenation 
at present, under the reasonable assumption that all minor 
mergers occur with gas-rich satellites. Thus, we expect that 
~ 10 — 30% of SO galaxies appear somehow rejuvenated by 
a minor merger at z = 0. 

On the other hand, we must keep in mind several caveats 
about this result. First, the minor-merger fraction esti¬ 
mated from cosmologic al simulatio ns is quite uncertain. 
Furthermore, IBertone fc Conselicel 12 0091 4 show that the 
major-merger fraction estimated from the same cosmolog¬ 
ical simulations is in partial disagreement with the data. 
The same issue might affect even the estimate of minor 
mergers. 

We can compare such estimates based on cosmolog¬ 
ical simulations with several observational signatures of 
minor mergers. Shell structures are one of the strongest 
scars leaved by merging episodes in ETGs. The litera¬ 
ture provides quite different estimates about the fraction 
of ETGs showing shells, since the shell pheno meno n likely 
depends on the galaxy environment (see e.g. lMarcum et all 
1200411 . The fraction of ETGs showing systems of shells 
range s from ~ 16% dMalin fc Carted 19831 Reduzzi et al.l 
Il996lf up to ~ 56% (according to Seitzer fc Schweizeil 1990. 
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who study shell frequency among 74 field ETGs). A high 
fraction of shells, 44% and 41% are found in ETG s by 
ISchweizer fe Ford! (1198511 and bv IColbert et al.l (|2001h . re¬ 
spectively. 

Recently, the ultra-deep survey MATLAS investi gated 92 
out of 260 ETGs in the ATLAS 3£I sample (|Puc et all 
!2014l l. Shell structures are detected in 21% of the ETGs ex¬ 
amine d. Furthermore, about half of the ETGs in lDuc et al.l 
: 201 11 sample show indications of minor merger (16%), ma¬ 
jor merger (12%), or interaction (22%). A large fraction 
of ETGs show perturbed/asymmetric halos (35%), streams 
(28%) and tails (17%). In summary, previous studies indi¬ 
cate that a significant fraction of nearby ETGs, including 
SO galaxies, may have suffered minor merger episodes. 

A possible estimate of the fraction of ETGs that un¬ 
derwent a recent SF episode may be obtained from mid 
infrared (MIR) spectra. Polycyclic aromatic hydrocarbons 
(PAH) visible at MIR wavelength s are c onsidered a clear 
signature of recent SF activity rtThielensl . 1 20081) . In ETGs, 
PAH s are detected with b oth normal and a nomalous ra¬ 
tios (iKaneda et al.l . I2008T) . IVega et all (|2010lf suggest that 
anomalous PAH ratios are due to the injection of carbona¬ 
ceous material in the inter-stellar medium by carbon stars, 
which are present in the stellar populations with ages in the 
range of a few Gyr dNanni et al.l . l2013h . MIR spectra should 
then be able to provide indications about both on-going 
and re cent (few Gyr) SF activity. Recently. iRanmazzo et all 
(j2013ll classified MIR Spitzer-IRS spectra in the nuclear re¬ 
gion (2 — 3 r e /8, where r e is the effective radius) of 99 ETGs 
according to galaxy activity. In general, elliptical galaxies 
are significantly more passive than SO galaxies. 46^{q% of 
elliptical galaxies and 20+7 1 % of SO galaxies show passively 
evolving spectra, i.e. spectra without an y signat ure of ei¬ 
ther emission lines or PAHs (see also iBressan et al.l . 120061) . 
PAHs are found in 47^®% of ETGs, but only 9lg% have 
PAH ratios typical of star forming galaxies. MIR spectra 
then suggest that a significant fraction of ETGs have been 
rejuvenated in their nuclear region in the last few Gyr. Our 
simulations show that both the ring and the nuclear region 
are interested by the rejuvenation of the stellar population. 

The aforementioned results are fairly consistent with the 
fact that about half of the SO galaxies in the local Universe 
show signs of rejuvenation, and suggest that a significant 
fraction of these episodes in nearby SO galaxies are induced 
by minor mergers with gas-rich satellites. 

Furthermore, lAnnibali et al.l (2007) modelled the 
nuclear Lick line-strength indices 65 n earby ETGs 
(|Ramoazzo et all 120051 : lAnnibali et all 120061 1 . By compar¬ 
ing the number of rejuvenated ETGs with the total number 
of galaxies in the sample, and by means of simple two- 
component single stellar population models, they estimate 
that the rejuvenation episodes do not involve more than 
25% of the total galaxy mass. This estimate, together with 
morphological signatures and MIR indication, witnesses 
that the role of minor mergers could be very important in 
the evolution of ETGs in general, and of SOs in particular. 

How efficient is gas transfer in minor mergers? Based on 
the WHISP survey (The Westerbork HI survey of irregular 
and spiral galaxies. lvan der Hulst et aLlteOPlU . ISancisi et alJ 
( 2008 ) estimate that ~ 25% of local galaxies show signs of 
possible minor mergers (which is in agreement with the 
aforementioned value from cosmological simulations), but 
find that the mass accretion rate is only ~ 1/10 of the 


SF rate. Recently. iDi Teodoro fc Fraternalil (1 2014h estimate 
that 22% of galaxies in a sub-sample of the WHISP have a 
dwarf companion, but the maximum gas accretion rate from 
these companions is about five times lower than the average 
SF rate of the sample (see also lHolwerda et akll201lD . These 
results suggest that minor mergers do not play a significant 
role in the total gas accretion bud get. On_the o ther hand, 
the galaxies consid ered by both ISancisi et al.l (|2008l l and 
IDi Teodoro fe Fraternalil (|2014l ) are gas-rich spiral galaxies. 

Spiral galaxies have their own reservoir of gas, and SF 
episodes might be triggered in several ways. In contrast, SO 
galaxies do not have a rich reservoir of gas of their own: for 
SO galaxies, the infall of a gas-rich companion might be the 
privileged way to refuel and re-activate SF. 

5. Conclusions 

In this paper, we investigated interactions between a SO 
galaxy and a gas-rich small satellite galaxy, by means of N- 
body/SPH simulations. The satellite galaxy is initially on a 
nearly parabolic orbit with respect to the primary galaxy. 
We ran six simulations, varying the impact parameter (10 
and 30 kpc), the inclination with respect to the plane of the 
SO disc (0, 45 and 90°), and the orientation of the orbital 
angular momentum of the secondary galaxy with respect to 
the spin of the primary galaxy (prograde and retrograde). 
In all simulations, the satellite galaxy is slowly disrupted 
by the interaction, its gas is stripped during each periapsis 
passage, and accretes onto the SO galaxy. 

We find that about 1/10 of the total gas mass of the 
satellite galaxy ends up in the central region of the SO 
galaxy. Warm dense gas rings form in most our simulations, 
indicating that minor mergers are a viable scenario for the 
formation of gas rings in SO galaxies and for the rejuvena¬ 
tion of SO galaxies. The radius of the ring is of ~ 6 —13 kpc, 
depending on the simulation, and the total mass confined 
in the ring is 1/100 — 1/10 of the initial gas mass of the 
satellite galaxy. 

We find that gas is stripped earlier in prograde runs than 
in retrograde ones, because tidal stripping is more efficient 
in prograde runs, where the relative velocity between the 
satellite galaxy and the particles in the SO disc is lower. 
Thus, prograde encounters build gas rings earlier, but these 
rings have a shorter life (1 — 3 Gyr) than those formed in 
either retrograde encounters or non-coplanar encounters. 
Rings formed in retrograde encounters and non-coplanar 
encounters can live unperturbed for more than ~ 6 Gyr. 

Furthermore, the rings keep memory of the orbit of the 
satellite galaxy: rings formed in prograde encounters are co¬ 
rotating with the disc of the SO galaxy, while rings formed in 
retrograde encounters are counter-rotating. Similarly, satel¬ 
lite galaxies whose orbit was coplanar with respect to the 
disc of the SO galaxy form co-planar rings, while satellite 
galaxies whose orbit was inclined by 45 and 90° with re¬ 
spect to the disc of the primary galaxy form rings that are 
inclined by 45 and 90°, respectively. 

This result has another crucial implication: rings form 
and grow even without bar resonances. In our simulations, 
the strongest and most long-lived ring is polar (run D): it 
forms with an inclination angle of 90°. 

Furthermore, the impact parameter appears to be a very 
important ingredient: if it is too large (30 kpc) and if the 
run is retrograde, the stripping is very inefficient and rings 
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do not form for the entire duration of the simulation (> 10 
Gyr). 

The formation of stellar shells surrounding the SO 
galaxy is also ubiquitous in our simulations, but such shells 
might be very faint: the mass density of stars in the shells 
is about three orders of magnitude lower than that in the 
disc of the SO galaxy. 

During the first periapsis passage, SF occurs mainly in 
the disc of the satellite galaxy, which is strongly perturbed 
by the interaction. As soon as the stripped gas reaches 
the centre of the SO galaxy, SF switches on even in the 
SO galaxy, while it is fast quenched in the satellite galaxy, 
whose initial gas was in large portion removed. The stripped 
gas powers an extremely long-lived (~ 8 Gyr) episode of 
SF in the central region of the SO galaxy. The SF rate is 
~ 10 -4 — 0.01 M 0 yr _1 , depending on the simulation. We 
expect that it can reach higher values for a larger initial 
gas mass of the satellite galaxy and/or if we account for 
the gas reservoir of the SO galaxy. SF and stellar feedback 
(supernova explosions) might play a strong role in the sup¬ 
pression of the gas rings: we find that gas rings disappear 
earlier in the runs where SF starts earlier and reaches the 
highest level (runs B and F). 

Finally, we discuss our results in view of the minor- 
merger fracti on deriv ed from cosmological simulations 
(jBertone fc Gonselicel . 12009 b from morphological signa- 
t ures (shel l s and streams) investigated in deep surveys 
(iDuc et all I2014T ), and from the presence of PAHs in the 
nuclear region of SO galaxies. We find that minor mergers 
might account for a significant fraction of rejuvenated SO 
galaxies. This result is a fundamental clue for our knowl¬ 
edge of rejuvenation processes in ETGs, but deserves fur¬ 
ther investigations, since it relies on a very uncertain esti¬ 
mate of the minor-merger rate in the local Universe, and 
on a limited sample of galaxies investigated with MIR and 
optical spectroscopy. 
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